The quantum Heisenberg antiferromagnet on the square lattice 
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The pure-quantum self-consistent harmonic approximation, a semiclassical method based on the 
path-integral formulation of quantum statistical mechanics, is applied to the study of the thermody- 
namic behaviour of the quantum Heisenberg antiferromagnet on the square lattice (QHAF). Results 
for various properties are obtained for different values of the spin and successfully compared with 
experimental data. 



We consider the quantum Heisenberg antiferromagnet 
on the square lattice (QHAF), whose Hamiltonian reads 



H = J J2 Si • ^ 



(1) 



<u> 



J is positive, the sum runs over all the couples < ij > of 
neighbouring sites on the square lattice, and the quantum 
operators Si obey the angular momentum commutation 
relations [S?, Sf] = iS^S i3 s a ^ with \Si\ 2 = S(S + 1). 

Several real compounds are well described, as far as 
their magnetic behaviour is concerned, by this model 
with S = 1/2 (La 2 Cu0 4 , Sr 2 Cu0 2 Cl 2 ), S = 1 (La 2 Ni0 4 , 
K 2 NiF 4 ) and S = 5/2 (KFeF 4 , Rb 2 MnF 4 ), and a con- 
sequently rich experimental analysis of the subject has 
been developed in the last ten years. From the theoreti- 
cal point of view, an equally rich reservoir of results, from 
both analytical and numerical approaches, is now avail- 
able; nevertheless, there are still many open questions, 
and different conclusions have been recently drawn by 
several authors [p]0||- 

To study the QHAF, we have used the pure-quantum 
self-consistent harmonic approximation (PQSCHA) 0, 
which is a semiclassical method based on the path- 
integral formulation of quantum statistical mechanics. 
Rs main feature is that of exactly describing the clas- 
sical behaviour and fully take into account the linear 
part of the quantum contribution to the thermodynam- 
ics of the system, so that the self-consistent harmonic 
approximation is only used to handle the pure-quantum 
nonlinear contribution. The fundamental goal of the 
PQSCHA is that of reducing the evaluation of quantum 
statistical averages to classical-like expressions involving 
properly renormalized functions, the fundamental one be- 
ing the effective Hamiltonian H c s- If j3 = T 1 , N is 
the number of lattice sites, Si is a classical vector on 
the unitary sphere (|si| = 1), and J d N s indicates the 
phase-space integral for a classical magnetic system, the 
quantum statistical average of a physical observable de- 
scribed by the quantum operator O turns out to be (O) = 
l/Z J d N s O exp(-0H eS ), where Z= J d N s exp(-/3W eff ) 
is the partition function. Both 7Y c ff and O depend on 
T and S, and the determination of their explicit form, 
starting from the expression of the original quantum op- 
erators, is indeed the core of the method Q . The effective 



Hamiltonian for the QHAF, i.e. relative to Eq. (JTJ), is 

= V 



JS 2 



= e 4 J2 + £(*) 



(2) 



where ,5 = 5+1/2,* = T/JS 2 and Q(t) is a uniform term 
that does not affect the evaluation of statistical averages. 
The renormalization coefficient 8 2 = 2 (t, S) < 1 is easily 
evaluated, for any given t and S, by self-consistently solv- 
ing two coupled equations ||. ^From Eq.(^) we see that 
the quantum effects leave the simmetry of the Hamilto- 
nian unchanged and introduce an energy scaling factor 
6 , naturally defining the effective classical temperature 



*eff 



t 



e\t,s) 



(3) 



that will enter all the PQSCHA results. The partition 
function, for instance, is Z = exp[— PG(t)]Z c \(t e g), where 
2 c \(t e f{) is the partition function of the classical model 
at a temperature i c g; O c i(£cff) will hereafter mean the 
value taken by the quantity O in the classical Heisenberg 
antiferromagnet at a temperature i c ff- 

The internal energy per site is easily found to be u(t) = 
9 4 (t, S)u c \(t e g) , while the correlation functions G(r) = 
(Si-Si+r), with i = (ij., ia) and r = (n, r 2 ) any vector on 
the square lattice, turn out to be 



G(r,i) = S 2 4 r G cl (r,t eS ) 



(4) 



the renormalization coefficients O 2 = 2 (t, S) are such 
that O 2 does not depend on r for large |r|, and O 2 = O 2 
for |r| = 1. From Eq.(|J), we find the PQSCHA expression 
for the staggered susceptibility x = X) r (~ ) ri+V2 G(r,t)/3 
to be 



S(S + 1) + S- 



E(- 



y 1+ r 2 04 Gci(rjtefi) 



(5) 



The PQSCHA result for the correlation length, de- 
fined by the asymptotic expression G(r) oc exp(— |r|/£) 
for large |r|, is £(t) = £ c i{t B s), meaning that the correla- 
tion length of the QHAF at a temperature t equals that 
of its classical counterpart at a temperature t e g . 

Once the problem has been reduced, by the PQSCHA, 
to a renormalized classical one, the ingredients needed to 



1 



obtain the temperature and spin dependent thermody- 
namic properties of the QHAF are the renormalization 
coefficients 9^(t, 5), whose evaluation is a simple matter 
of a fraction of second on a standard PC, and the tem- 
perature dependence of the corresponding properties of 
the classical model, typically obtained by classical Monte 
Carlo simulations 0. 

In the following we will focus our attention on the 
staggered susceptibility and the correlation length, as 
experimental data for these quantities are available for 
various compounds. Such compounds are usually char- 
acterized by a crystal structure in which the magnetic 
ions form parallel planes and mainly interact if belong- 
ing to the same plane; a weak interplanar interaction is 
responsible for a low-temperature 3D transition, and it 
introduces also an anisotropy term. Keimer et al. |?]] 
have shown that in the classical limit, and to one- loop 
level, the relation between £ in presence of the anisotropy 
term, and £o of the fully isotropic model, is given by 
£ = £o/(l — o^q) 1 / 2 , where a is a parameter describing 
the relative strength of anisotropy; following Lee et al. || 
we shall employ the above formula (with some refine- 
ments H which lead to substitute a with its renormal- 
ized counterpart a c s ) to compare our PQSCHA results 
with experimental data. 

In Figs. [I] and || we present our results for the correla- 
tion length of the QHAF for 5 = 5/2, and for the same 
quantity and the staggered susceptibility for S = 1. For 
5 = 5/2 the experimental data refer to Rb2MnF4 || and 
KFeF4 |)|: the anisotropy term is seen to be very well 
described by the approach described above. For 5 = 1 
QMC |]l0| and experimental data for the two compounds 
La 2 Ni0 4 Q and K 2 NiF 4 Q are reported. 

We may thus confidently conclude that the thermody- 
namic behaviour of the QHAF is properly described by 
the PQSCHA, i.e. in term of a renormalized classical 
Heisenberg antiferromagnet. The easy-axis anisotropy, 
sometimes crucial to analyse the low-temperature exper- 
imental data, has been considered following Keimer et 
al. 0, in the PQSCHA framework. Finally, we would 
like to recall that our results help to shed some light on 
the reasons of the failure of the theory based on the non- 
linear a model approach []l3|, in describing the QHAF 
for 5 > 1 in the temperature region where experimental 
data are available. 
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FIG. 1. Correlation length £ versus t for S = 5/2: Lines 
are the PQSCHA results for the QHAF with the bare a (dot- 
ted), the renormalized a e ff (dashed) and without anisotropy 
(full). Experimental data for Rb2MnF4 || (squares) and 
KFeF 4 § (circles). 
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FIG. 2. Correlation length £ and staggered susceptibility 
3x* = 3Y/S 1 versus t for S = 1: Experimental data for 
La 2 Ni0 4 @ (squares hand K 2 NiF 4 [[2) (circles) ; QMC data 
(triangles) from Ref. |UJ. 
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